On the energy-minimizing steady states 
of a thin film equation 



Almut Burchard 
O ■ University of Toronto 

almut @math . utoronto . ca 

Oh 1 

go ; 

Marina Chugunova 

■ University of Toronto 

chuqunom@math . utoronto . ca 

< 

^ ■ Benjamin K. Stephens 

University of Washington 

benstph@math . Washington . edu 

> 

On: March 9, 2010 

O 

ON 

o 
o 

Abstract 

Steady states of the thin film equation ut + [u 3 {uqqq + a 2 ug — sin 9)]q = are considered 
on the periodic domain Q = (—ir,ir). The equation defines a generalized gradient flow for an 
energy functional that controls the H 1 (£l)-norm. The main result establishes that there exists for 
each given mass a unique nonnegative function of minimal energy This minimizer is symmetric 
decreasing about 6 = 0. For a < 1 there is a critical value for the mass at which the minimizer has 
a touchdown zero. If the mass exceeds this value, the minimizer is strictly positive. Otherwise, it 
is supported on a proper subinterval of Q and meets the dry region at zero contact angle. A second 
result explores the relation between strict positivity and exponential convergence for steady states. 
It is shown that positive minimizers are locally exponentially attractive, while the distance from a 
steady state with a dry region cannot decay faster than a power law. 
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1 Introduction 



Degenerate fourth order parabolic equations of the form 

Ut + V • (/VA«) + lower order terms = 

are commonly used to model the evolution of thin liquid films on the surface of a solid. Here, u(x, t) 
describes the thickness of the fluid at time t at the point x, the fourth derivative term models surface 
tension, and the exponent n > is determined by the boundary condition between the liquid and the 
surface of the cylinder. A particularly interesting case is n = 3, which corresponds to a "no-slip" 
boundary condition. 

In this paper, we study the equation 

U t + [u n (uggg + G? Ug — SU1 6) + Uu\ Q = , 8 G f2 = ( — 7T, 7r) (1.1) 

with periodic boundary conditions. For n = 3 and a — 1, this describes the evolution of a thin liquid 
film on the outside of a horizontal cylinder that rotates slowly about its axis, see Figure 1 . The film is 
assumed to be uniform along the axis of the cylinder, and its thickness at time t and angle 6 (measured 
from the bottom) is given by the function u(9, t). In Eq. (11.11) . the first summand in the parentheses 
models surface tension, the next term is a correction due to the curvature of the cylinder, and the third 
term models gravitational drainage. The last term models convection due to rotation. We have scaled 
the units of length and time so that the surface tension and gravitational terms appear with coefficient 
one; the coefficient a > is a geometric constant, and u is proportional to the speed of rotation. Here, 
we will study the the non-rotating cylinder with u = 0. We are mainly interested in the case where 
n = 3 and a—1, but find it illuminating to also consider other values of n and a. 
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Figure 1 : Liquid film on the outer surface of a rotating horizontal cylinder in the presence of gravity. 

The model in Eq. (11.11) (with n = 3 and a = 1) was first derived and studied by Pukhnachov 
[|24ll25l . Note that Pukhnachov uses slightly different coordinates (with 8 = ^ at the bottom of the 
cylinder), and a different scaling (where u — 1). The derivation relies on the lubrication approxima- 
tion; it assumes that the liquid film is very thin compared to the radius of the cylinder, and that the 
rotation is slow enough to neglect centrifugal forces. Pukhnachov 's model refines an earlier model of 
Moffatt [12711 by including surface tension. The second order curvature term in the equation is rem- 
iniscent of a porous-medium equation, but appears with the opposite sign, resulting in a long-wave 
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instability J6]|. Interesting numerical and asymptotical analysis along with numerous open questions 
can be found in HIHEI. 

A fundamental technical problem in thin film equations is well-posedness, i.e., to show that non- 
negative initial data give rise to unique nonnegative solutions that depend continuously on the data. 
The difficulty is that solutions of fourth order parabolic equations generally do not satisfy maximum 
or comparison principles, and linearization leads to semigroups that do not preserve positivity. To 
give a simple example, the function u(x, t) — 1 + tcos(x), which solves Eq. (11.11) in the linear case 
n = with a = 1 and u — 0, ceases to be positive for t > 1. Thus, positivity of solutions of thin 
film equations is a nonlinear phenomenon. Other relevant problems in the area concern long-term 
behavior, finite-time blow-up, and the interface between wet (u > 0) and dry [u = 0) regions. These 
problems have been studied rigorously in a vast body of papers since the pioneering article of Bernis 
and Friedman 01, see for example [J3l El [fE [TOl |22]| and references therein. There is an even larger 
literature that studies the properties of physically relevant solutions through asymptotic expansions, 
numerical analysis, and laboratory experiments. 

Bernis and Friedman proved that initial- value problems for thin film equations can be solved in 
suitable classes of nonnegative weak solutions. An important technical contribution was their use of 
an entropy functional that decreases with time along solutions. A few years later, the so-called a- 
entropies were independently discovered by Bertozzi and Pugh [5] and by Beretta, Bertsch, and dal 
Passo [3]. Since then, other families of entropies have been found jH[T3l[T6l|. Entropy functionals are 
the basis for results on short- and long-term existence, positivity, finite speed of propagation, regularity, 
blow-up, and the long-time behavior of solutions. 

For < n < 1, well-posedness and convergence to steady states have recently been established by 
treating Eq. (11.11) as a gradient flow on a space of measures endowed with the Wasserstein distance [20, 
1221 . where the exponent n appears as a mobility parameter [9]. These gradient flow techniques also 
take advantage of energy and entropy dissipation. However, Wasserstein distances with mobility n > 1 
are not well understood, and well-posedness remains an open problem. 

By a solution of Eq. (11.11) we mean a nonnegative function u E L 2 ((0, T), H 2 (Vl)) that satisfies 



for every smooth test function with compact support in Q x (0, T). This agrees with the strongest 
notion of generalized solutions from [HUH. For a class of equations that includes Eq. (|l.ll) with n = 3, 
long-time existence of such solutions was recently proved in [10]. These solutions are widely believed 
to be unique. 

Questions about steady states appear in many applications. When do steady states exist, when 
are they uniquely determined by their mass, are they strictly positive or do they exhibit dry regions, 
and under what conditions are they stable? Do steady states attract all bounded solutions? When can 
we expect exponential convergence? Linearizations of Eq. (11.11) about steady states were examined 
analytically and numerically in lfTT][T2ll . 

For u = 0, Eq. (11.11) defines a generalized gradient flow for the energy 



in the sense that u t = \u n (ff )J e - Here, ^ denotes the L 2 -gradient of E. This implies the dissipation 





(1.2) 
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estimate 



d_ 

dt 



*«■•'» = -/„«"(£) 



dO < 0. 



(1.3) 



The subject of this paper are the minimizers of E on the set of nonnegative 2n -periodic functions of a 
given mass 



and their role in the dynamics. Note that for a = 1, E is convex but not strictly convex. For a < 1, 
the functional is strictly convex, and has a unique critical point on Cm, which is a global minimizer. 
For a > 1, it is not convex, and we may expect multiple critical points. 

We will show that E has a unique minimizer on Cm for each value of a and each mass M > 0, see 
Theorem [T] These minimizers may have dry regions; in that case, the contact angle of the fluid film is 
zero, see the Figure 2. In particular, our result establishes the existence of zero contact angle steady 
states for Eq (|l.ll) if M (1 — a 2 ) < 2tt. Our proof relies on symmetric decreasing rearrangements and 
the first variation of the energy. 

The minimizers are time-independent solutions of Eq. (Il.lt with uj = 0. Additional steady states 
may arise for a > 1 as saddle points of the energy. For any value of a and M, there is also a continuum 
of steady states with have non-zero contact angles, analogous to the steady states in 031, whose role 
in the evolution remains unclear. 

We expect that for a < 1, the unique energy minimizer should attract all solutions on Cm as 
t — >■ oo. Unfortunately, in the absence of a proper well-posedness theory, Lyapunov's theorem is not 
sufficient to support this expectation. In Theorem |2] we provide partial results in that direction. If 
the energy-minimizing steady state is strictly positive, then it exponentially attracts all solutions in a 
neighborhood. On the other hand, for n > |, a steady state that has a dry interval of positive length 
cannot be exponentially attractive. In particular for n = 3 and a = 1, the distance between the solution 
and the minimizer decays no faster than £~3. Our proof combines energy and entropy inequalities in 
the spirit of E [51 HIH. 

All our results are easily adapted to the long-wave stable case where the sign of the second order 
term is reversed. In that case, the energy-minimizing steady state is strictly positive and locally ex- 
ponentially attractive so long as M (1 + a 2 ) > 2ty. For M(l + a 2 ) < 2ir, the energy minimizer has 
a dry interval of positive length, which it meets at zero contact angles. This contrasts with a theorem 
of Laugesen and Pugh that excludes zero contact angle steady states for the corresponding thin film 
equation without the sine term [fT51 . 

2 Identification of energy minimizers 

We begin by showing that for every M > there exists a function u with mass M that minimizes the 
energy. The first lemma provides the necessary global bounds on the functional. 

Lemma 1 (Lower bound on the energy.) E is bounded from below and coercive on C M . 

Proof. Using that u is nonnegative and has mean we estimate 





u\\oo < 7T + V^IKIU 2 • 
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It follows that 



E(u) > \ (\\u\U - ^(1 + a 2 )]' - ^a 2 (2 + a 2 ) - M. (2.1) 



which is clearly bounded below. A similar estimate shows that E grows quadratically as | \u\ \ H \ — > oo. 

□ 

The lemma implies that minimizing sequences are bounded in H l . Passing to a subsequence, we 
can construct a minimizing sequence {uj}j>i that converges weakly in H 1 and strongly in L 2 to some 
function u in Cm- Since E is weakly lower semicontinuous on H 1 , u is the desired minimizer. We 
next describe some properties of the minimizers. 

Lemma 2 (Symmetry.) Minimizers of E on Cm are symmetric decreasing about 9 = 0. 

Proof. For u G Cm, let -u # be the unique symmetric decreasing function of 9 that is equimeasurable 
to u. Classical results about symmetric decreasing rearrangement ensure that u* 6 Cm, and that 

| |it#| \ L i = \\u\\i,2 , | \uf\ \i? < \\uo\\ip. , / u* cos 9 d6 > / ucosOdO, (2.2) 

Jn Jq 

see |[23]|. It follows that 

E{u*) < E{u) . 

If u is a minimizer, then E(u#) = E{u), and in particular, the third inequality in Eq. (12.21) must 
hold with equality. Since the cosine is strictly symmetric decreasing, this forces u to be symmetric 

decreasing as well [fTTl Theorem 3.41. □ 

The Euler-Lagrange equation for the minimizer is given by 

u ee + a 2 u + cos 9 = X on {9 e Q | u(9) > 0} , (2.3) 

where A is a Lagrange multiplier associated with the mass constraint. By considering the first variation 
of E with respect to nonnegative functions that need not vanish outside the support of u, we see that 

uqq + a 2 u + cos 9 < A on Q 

as a distribution. This suggests that the first derivative of a minimizer should vanish at the boundary 
of its support, i.e., the film meets the surface of the cylinder at zero contact angle. The next lemma 
confirms this suspicion. 

Lemma 3 (Zero contact angle.) Let u be a minimizer of E on Cm- If u has its first positive zero at 
9 = r, then uoij^) = and ueeij ) > 0. In particular, u G C l,1 {Vl). 

Proof. By Eq. (12.31) . u has one-sided derivatives of arbitrary order at r, and uq{tJ) < 0. Suppose that 
M 6 i(r_) < 0. We will modify u to construct a new valid competitor with lower energy. If r < it, set 
v (9) = u((fi(9)), where (p : Vt — > is the bi-Lipschitz map defined by 

r-e < 9 < r + e, 
t + s<9<t + 2s, 
otherwise . 
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Figure 2: On the left: Steady states for a = 1 and initial data uo = 0.5, 1,2. As the mass is getting smaller, the 
minimizer becomes more concentrated, and as mass goes to infinity the support tends to [—it, tt}. On the right: 
Energy-decreasing variations resulting from non-zero contact angles. 

Choose e > small enough so that u vanishes on [r, r + 2e\. Then v vanishes on [r + e, r + 2e\. The 
difference between the leading terms in the energy integrals is given by 

\j/ode- l -f^ e d0=\ jT* vl 4 de = -|«, (r_) a + 0(e 2 ) , 

see Figure 2, top right. Since the remaining terms contribute only corrections of order 0(e) 2 to the 
energy difference, it follows that 

E(v)-E(u) = - £ -u (r^) 2 + O(e) 2 . 

If r = Ti, the same estimate holds (by the symmetry of u) for v(9) = u(min{|0|, n — \s}), see 
Figure 2, bottom right. In either case, the mass of v is M' = J n v d6 = M + 0(e) 2 . We finally 
set w(9) = jpv(9), which has the correct mass and satisfies E(w) = E(v) + 0(e 2 ) < E(u) for e 
sufficiently small, a contradiction. We conclude that u e (r_) = 0, proving the first claim. 

To prove the second claim, we analyze the sign of the first non-vanishing left derivative of u at 
t. Clearly, uee(r^) > 0. Suppose that uee(r-) = 0. Differentiating Eq. (12.31) . we obtain u eee (T-) = 
sinr > 0. Since r is the first positive zero of u, the derivative u from the left cannot be positive, and 
so 9 = 7r is the only possibility. Differentiating once more, we obtain u edee (ir_) = — 1, which is the 
wrong sign for u to have a minimum at n. It follows that u eg (r_) > 0, as claimed. 

□ 

We are now ready to compute the minimizers of E on Cm explicitly. A particular solution of the 
Euler-Lagrange equation in Eq. (12.31 ) is given by 

{— ~9 sin 9 , a = 1 , 

^(cos^-l^cosM)), a^l, (2A) 

and the general solution can be represented as 

u(9) = A cos(a9) + B sm(a9) + ^ + u°(9). (2.5) 
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The following theorem summarizes our results. The statement is illustrated in Figure 3. 




Figure 3: On the left: Numerical steady states for a = 0.5, 1, 2 with initial data uq = 3. On the right: Mass 
versus half length of the compact support for a = 2,1, 0.5. 



Theorem 1 (Description of the energy minimizers.) Let E be the energy functional in Eq. ([772]). For 

each M > 0, E has a unique nonnegative minimizer of mass M. The minimizer is strictly symmetric 
decreasing on its support. It is of class C 1 ' 1 and depends continuously on M in C 1 ' 1 . It increases with 
M in the sense that for any pair of minimizers Ui, u 2 of mass Mi, M 2 , 

Mi < M 2 =► ui(9) <u 2 (9), 9 e support (m) . 

• If M{\ — a 2 ) > 2ir, the minimizer is strictly positive and given by 

M 1 

u{9) — - — H -cos 9; (2.6) 

Z7T 1 — cr 

• if M{1 — a 2 ) < 2n, the minimizer is given by 

u {9) = A(r)(cos(a6>) - cos(or)) + u°(9) -u°(r) , |0| < r (2.7) 

for some r with max{«, 1} r < 7r, and vanishes for \9\ > r. The coefficient is determined by 
A(t) = a ^(aT) ' wnere u ° i s tne special solution from Eq. (\2.4\) ; 

• ifM{\ - a 2 ) = 2tt, then Eq. $273 for r = rc coincides with Eq. ri26l) . and u(9) = 

Proof. Fix M > 0. By Lemma [T} there exists a minimizer u of mass M, and by Lemma [2l it is 
symmetric decreasing about 9 = 0. If the positivity constraint is not active, then Eq. (12.51) holds for all 
9 E £1. Since the minimizer is smooth, periodic, and has mass M, we conclude that a < 1 and Eq. (12.61) 
holds. In order for u to be nonnegative and symmetric decreasing we must have M{1 — a 2 ) > 2n. In 
that region, u is clearly strictly increasing in 9. 

If, on the other hand, the positivity constraint is active, then the minimizer u is positive on some 
interval (— r, r) and vanishes for \9\ > r. By Lemma[3l u E C 1 ' 1 ^) and ite(±r) = 0. On its support, u 
is given by Eq. (12.51) . Since u and u° are even, B = 0. The Dirichlet condition at r allows to eliminate 
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A, the Neumann condition determines A, and we find that Eq. (12.71) holds. If max{a, 1} r < it, we 
claim that u is indeed nonnegative, symmetric decreasing in 9, and strictly increasing with r. To see 
this, we differentiate Eq. (|2.7I) . and use Lemma[3]to obtain 

dA . . . , 

— ■ a sin ar = —ug T {T; t) = uge{r; r) > . 

By the chain rule, and using once more that ug(r; r) = 0, we have 

dA 

u T (9; t) = — ■ (cos(a6>) — cos(ar)) > for \9\ <r . 
dr 

Since u vanishes identically when M = 0, this confirms that u is positive and strictly symmetric 
decreasing for \9\ < r. We use that uq(t] t) = to compute 

dM dA [ T 

— — = — — / cos(cw) — cos(ar) d9 > , 
dr dr J_ T 

and infer that we can solve for r = r(M) as a strictly increasing smooth function of M. By the chain 
rule and the inverse function theorem, 

d .„ , , cos(a60 — cos(ar) 

— u(9; t M = > , 

dM J_^_ cos(afc' / ) — cos(ar) ad' 

establishing the desired continuity and monotonicity of u with respect to M in the range where 
Eq. (12.71) is valid and max{a, 1} r < it. 

We need to determine the ranges where Eq. (12.61) and (12.71) hold. For a < 1, E is strictly convex. If 



M > yz^2> the function defined by Eq. (|2.6I) is nonnegative and provides the unique minimizer of E 
on Cm- If M < jz^, the positivity constraint is active, and we compute from Eq. (12.71) that M — > 
as r -> and M — > as r — >• 7r_. Continuous dependence on M follows, since Eq. (12.71) agrees 
with Eq (fTll) at M = 

For a = 1, E is convex, but not strictly convex on Cm- The positivity constraint is active, because 
E is not bounded below without it; for instance, E + tcos9) = —nt. Note that we must have 
r < 7r, because the particular solution u°(9) = —^9sm9 from Eq. (|2.4I) cannot be continued as a 
differentiable periodic function across 9 = n, in violation of Lemma [3] It is easy to check from 
Eq. (12/71) that M -)■ as r -)> 0, and M ->■ oo as r -)> 7r. 

For a > 1, the energy is a non-convex quadratic function on Cm, and hence the positivity constraint 
is active and u is given by Eq. (12.51) on some interval (— r, r). Upon closer inspection of Eq. (|2.5I) . we 
see that ar < 1, since otherwise u fails to be symmetric decreasing. Since M — > as r — > and 
M — )■ oo as r — > a" 1 ^, the theorem follows. □ 



3 Convergence to minimizers 

In this section, we will prove a lower bound on the speed at which solutions of Eq. (11.11) with u = 
can converge to critical points on the boundary of the positive cone. We establish this bound for two 
classes of solutions: For strictly positive, classical solutions when n > |, and for the strong generalized 
solutions constructed for n = 3 in [10, Section 3]. 



8 



Our bound uses the entropy method, applied to the functional 



S(u) = I u-P d6 , 
'n 



(3.1) 



where /3 = n — |. Strictly speaking, S 1 is not an entropy for Eq. (11.11) , because it may increase as well 
as decrease along solutions. One of the reasons is that the porous medium term a 2 (u n u e ) e appears 
in Eq. with the unfavorable sign. Still, the standard entropy methods yields a useful differential 
inequality for S. 

Lemma 4 (Entropy inequality for classical solutions) Fix n > § and let S be given by Eq. ri3.il) 



with (3 = n 



For every strictly positive classical classical solution u of Eq. di.il) , there exist 



constants S and K such that 



S(u(;t))<S + K t 



(3.2) 



Proof. Let So = S(u(-, 0)) and E = E(u(-, 0)) be the initial values of the entropy and energy, and set 
c n = (n — |)- We will show that 



d_ 



S{u{;t))<K , 



where 



K = c n 



Ma 2 
4 



M M 2N 



2(E + M) M 5 



7T 



4tt 2 



a 2 (2 + a 2 



2^ 



7T 



The claim then follows by integrating along the solution. 

To see the differential inequality, we use the Eq. (11.11) and integrate by parts, 



u 2 uqUqqq d6 + a u 2 u 9 d6 — u 2 ug sin 6 d6 . 
n Jn Jn 



n dt 

The first summand in Eq. (13.31) we integrate again by parts, 



u zueUeeedO = - 
n Jn 



-I 2 7/1 1 / -3 2 1 

it 2-u dU + ■ ' 



u 2 u e u e e d6 =: —A + -B 



(3.3) 



and integrate by parts once more to see that B = \ J u zujdO. After collecting terms, the first 
summand becomes 



-A + B--B 
2 



u a ^u ee - 1 u 2 ^j d6 = -4 J u* ( (u2) ee ) 



The second summand in Eq. (|3.3I) we also integrate by parts, 



a 2 u zu 9 d6 = 2a / u(u^) ee d9, 
>n Jn 
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and the third summand we rewrite as 



u 2 ug sin 6 dO = 2 / u 2 cos 8 d6 . 
n Jn 



Inserting these identities into Eq. (13.31) and completing the square we arrive at 



dS(u) I f 1 / 1 a 2 i\ 2 a 4 



. ix* I 2(«3) M - _ u s d6 + — u*d8 + 2 utcosBdB} (3.4) 
[ Jn \ 2 / 4 y c 

< c n |^M v /||n|| 00 + 2v / M^| . 

The claim follows from Eq. (12.11) of Lemma [H □ 

The entropy inequality in Eq. (13.21) holds also for many types of weak solutions of Eq. (11.11) that 
are obtained as limits of classical solutions of suitable regularizations. We demonstrate this for the 
strong generalized solutions in the case n = 3 that were recently constructed by Chugunova, Pugh, 
and Taranets in IfTOl Theorem 2]. For e > 0, consider the regularized equation 

u t + [f e {u) (u eee + a 2 u e - sin 9)] e = , 8eQ, (3.5) 

where f £ {z) = j^r e - This equation has strictly positive classical solutions positive initial data in 
if 1 (f2) (see IfTOl Lemma 3.4]), and the energy in Eq. (11.21) is dissipated. 

The key step is to extend the entropy method to the regularized equation. Set 

SJu) = / s £ (u)dx, (3.6) 



where s e {z) = z 2 (1 + ~ez 1 ) is chosen so that s"(z)f £ (z) = c 3 z 2 . Then ^ satisfies the entropy 
identity in Eq. (13.41) along solutions of Eq. (13.51 ), and hence 

S E (u(;t))-S s (u(;0))<K t, 

with the same constant as in LemmalU As e — > along a suitable subsequence, solutions of Eq. (13.51 ) 
converge uniformly to solutions of the original problem in Eq. (11.11) . and the values of the energy and 
entropy also converge for all times, provided they are finite at t = 0. Thus the entropy inequality 
in Eq. (|3.2I) remains valid in the limit. This can be used to establish additional regularity properties: 
Using a suitable subsequence where ugg converges weakly in L 2 {Vt x (0, T)) one can show that for 
every T > 0, 

J j [2{vh)g B - y u ^) d8dt< S{u(-,0)) + K o T < 00, 

and conclude that j Q ((u^)gg^ d9 is finite for almost every t > 0. 

Our final result concerns the dynamics near the energy minimizer, see Figure 4. 

Theorem 2 (Bounds on the speed of convergence.) Consider Eq. ( Ii.il) with u = 0, and set f3 = 

n — |. Let u be a solution of mass M that dissipates energy and satisfies Eq. ( 13.21) with constants So 
and Kq, and let u* be the energy -minimizing steady state of the same mass. 
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Figure 4: Evolution of a solution with a = 1, n = 3 and initial data no = 1. On the left: Time shots of 
the numerical solution at t = 0, 10~ 2 , 10 _1 , 1, 10, 10 2 , 10 3 . On the right: L 2 -distance of the solution from the 
energy minimizer. The dashed line shows the bound from Theorem [2 



• Ifu* is strictly positive, set p,— (1 — a 2 ) (min u*) n and Eq = - % — (min u*) 2 . For each solution 
with E(u(-, 0)) < E(u*) + Sq there exists a constant K\ such that 

\H-,t) -u*\\ H i < ffie - "*; 

• ifu* vanishes on an interval of positive length L, and n > |, then 

\\u(-,t)-u*\\ 2 >L 1+ ^(S + K t)-^; 

• ifu* vanishes quadratically at a point, and n > |, then there exist positive constants K 2 and K 3 
such that 

|K-,t)-W*|| 2 > (#2+^3*)"^- 

Proof. If u* is strictly positive, then a < 1 by Theorem [Q and E is strictly convex. The Taylor 
expansion of E about u* terminates after the quadratic term, because E itself is quadratic, and the 
linear term vanishes because u* is a critical point where the positivity constraint is not active, and 
therefore 

E{u) = E{u*) + i I (u- u*) 2 d - a 2 {u - u*) 2 d6 . (3.7) 
2 Jn 

Since u and u* have the same mass, the energy difference dominates the Z/ 1 -distance, 

E(u) — E(u*) > 1 -^\\(u-u*) e \\l > (1 ~ 2 Q2)7r lk-^ll 2 oo- 
Set e = E(u(-, 0)) — E{u*) < e , then min u > (1 — — ) mmu* > 0. This means that u is a strictly 

■2 



positive, classical solution that can be differentiated as often as necessary. We compute the L -gradient 

6E 
Su 



of E from Eq. (13.71) as ^ = —(u — u*)ee — ct 2 (u — u*). By Eq. (11.31) the energy is dissipated at rate 



^-E(u(;t)) = - / ' u n [{u-u*) ee + a 2 {u-u*)\ 2 e de 
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< -(rninw)™- / [(u - u*) eee + a 2 (u - u*) e ] 2 d6 

Jq 

< -2 ^(l-i-)minu*^ (1 - a 2 ) (E(u) - E(u*)) . 

In the last step, we have used Parseval's identity to rewrite the integral and the energy difference in 
terms of the Fourier coefficients of u — u*, and estimated the Fourier multipliers by 

p 2 (p 2 -a 2 ) 2 >(l-a 2 )(p 2 -a 2 ), (p ? 0) . 

Exponential convergence of the energy follows from Gronwall's lemma. Since min u converges expo- 
nentially to min-u*, we conclude that E(u(-, £)) — E{u*) < Ke~ 2ilt for some constant K, and the first 
claim follows. 

If u* vanishes on an interval of length L > 0, we apply Jensen's inequality to the convex function 
y i — y y~ 2 on this interval to obtain 

S{u{-,t)) > L 1+ ^ \\u{-,t) -u*\\2 P . 

The second claim follows by using the bound on the entropy in Eq. (13.21) and solving for the distance 

IK-,*) - M*|| 3 . 

If u* vanishes quadratically at a point, we consider the interval of length L centered at that point 
and obtain with the same calculation as for the second case that 

\\u(-,t) - U*\\ 2 > ||lf(-,t)/|0|> T ||2 - ||w*/|0|>r||2 

> L^(S + K ty^ -O(Lf). 

The proof is completed by choosing L = e(S + K t)~~ for e > sufficiently small. □ 

Consider specifically the case of Eq. (11.11) where n = 3, a = 1, and to = 0. By Theorem [H the 
energy minimizer vanishes on an interval of length L = 2(n — r) > 0, and by Theorem |2l solutions 

2 

cannot converge to the minimizer more quickly than t~ 3 . We suspect that the distance from the steady 
state actually decays with t~3. 

This conjecture is supported by simulations, and by analogy with aggregation processes where 
the convergence to states is governed by power laws. In such processes, the speed of convergence 
is limited by the rate at which mass can be transferred from a region of low density to the region of 
accumulation. One example is the Lifshitz-Slyozov cubic law, which describes late-stage grain growth 
in alloys and the evaporation-condensation mechanism in supersaturated solutions lfT8l . A second 
example is the separation of a water drop from a tap, shown on the left of Figure 5. The drop slowly 
draws water from the tap through a thin neck. The right side of Figure 5 shows the accumulation of 
the mass of a solution of Eq. ( 11.11) at the bottom of the cylinder. The large drop at the bottom grows by 
slowly pulling mass from the small drop on the top through the thin bridges that connect them. 
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Figure 5: On the left: 
by Elliott Minns lfl9l . 
parameters n = 3, a ■ 



:: Separation of a water drop at a tap. Photograph "Perfect Back and White — Water drop", 
. On the right: Matlab simulation of the long-time behavior of a solution of Eq. (11.11) with 
= 1, and uj = 0. The thickness of the film is exaggerated to emphasize the shape. 
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